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An  Assessment  of  Potential  Nonlinear 
Circuit  Models  for  the  Characterization 
of  Resonant  Tunneling  Diodes 
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Abstract- The  intrinsically  fast  process  of  resonant  tunneling 
through  double  barrier  heterostructures  along  with  the  existence 
of  negative  differential  resistance  in  the  current-voltage  char¬ 
acteristic  of  these  structures  has  led  to  their  implementation 
as  sources  for  high  frequency  electromagnetic  energy.  While 
sources  based  upon  resonant  tunneling  diodes  (RTD  s)  have 
produced  frequency  of  oscillations  up  to  712  GHz,  only  microwatt 
levels  of  performance  has  been  achieved  above  100  GHz,  Since 
stabilitv  criteria  plav  a  critical  role  in  determining  the  deliverable 
power  of  any  oscillator,  a  physically  accurate  equivalent-circuit 
model  for  the  RTD  is  extremely  important  for  optimizing  the 
dynamics  of  the  device-cavity  package.  This  study  identifies  a 
distinct] v  new  equivalent  circuit  model  for  characterizing  the 
modes  of  oscillation  in  RTD-based  sources.  Specifically,  in  order 
to  exhibit  the  fundamental  self-oscillations  and  the  overall  I- 
V  characteristics  (plateau  structure  and  hysteresis)  observed 
experimentally,  an  accurate  circuit  model  of  the  RTD  must 
incorporate;  (i)  a  quantum-well  inductance  which  directly  chokes 
the  nonlinear  conductance  and,  (ii)  a  nonlinear  access  resistance, 
associated  with  the  accumulation  of  charge  in  the  injection  region 
of  the  double  barriers,  with  a  nonlocal  dependence  on  the  bias 
across  the  double  barrier  structure. 


I.  Introduction 

APPLICATION  of  resonant  tunneling  diodes  (RTD)  to 
high-frequency/high-speed  electronics  is  attractive  due  to 
the  fundamental  physical  principles  upon  which  they  operate. 
The  RTD,  as  was  first  demonstrated  by  Chang  et  al.  [1], 
exhibits  the  purely  quantum  mechanical  process  of  tunneling 
where  electrons  transverse  the  double-barrier  structure  via 
the  quasibound  energy  state  of  the  well.  Since  this  tunneling 
occurs  through  the  quantum-well  state(s)  which  first  becomes 
aligned  and  then  passes  by  emitter  states  in  the  source  [2] 
during  a  bias  sweep,  the  RTD  yields  a  negative  differential 
resistance  (NDR)  current  characteristic.  Furthermore,  since 
tunneling  is  an  intrinsically  fast  process  the  RTD  current  can 
respond  very  rapidly  to  changes  in  applied  bias.  In  fact,  RTD’s 
have  been  shown  experimentally  to  have  detection  capabilities 
up  to  2.5  THz  [3]  and  have  been  implemented  as  oscillators  up 
to  712  GHz  [4].  Therefore,  the  remarkable  fast-acting  NDR  of 
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RTD’s  make  them  excellent  devices  for  very  high-frequency 
electronic  applications. 

RTD’s  have  been  extensively  investigated  as  sources  of 
terahertz  (1012  Hz)  electromagnetic  energy.  Unfortunately, 
RTD-based  power  sources  have  only  demonstrated  microwatt 
levels  (i.e.  <50  iiW)  of  performance  above  100  GHz  [5]. 
Thus,  RTD’s  have  been  outperformed  in  this  frequency  band 
by  more  traditional  sources  (e.g.  IMPATT’s,  TUNNET’s,  etc.) 
and  the  power  levels  achieved  are  not  sufficient  for  many 
practical  applications  [6].  The  high  frequency  potential  of 
RTD’s  has  concentrated  much  effort  toward  the  theoretical 
optimization  of  these  structures  using  a  number  of  material 
systems  [7].  While  quantum  simulations  [8]-[10]  are  necessary 
to  derive  the  fundamental  characteristics  of  isolated  structures, 
equivalent  circuit  models  are  usually  employed  to  determine 
the  power  and  stability  limitations  of  an  RTD  within  an  actual 
resonating  cavity  or  circuit.  To  date,  such  estimations  of  RTD 
performance  have  been  based  upon  linear  oscillation  theory 
and/or  have  utilized  simple  small-signal  diode  models  [11]- 
Obviously,  the  results  hinge  upon  the  equivalent  circuit  model 
used  in  the  analysis.  Hence,  efforts  to  optimize  the  power 
performance  of  RTD-based  sources  can  be  curtailed  if  the 
choice  and/or  approach  to  circuit  modeling  is  not  an  accurate 
representation  of  the  dynamics  of  the  device-cavity  package. 
Therefore,  the  identification  of  a  equivalent  circuit  model 
which  closely  matches  the  physical  dynamics  is  extremely 
important  in  the  design  and  optimization  of  RTD  oscillators. 

A  phenomenon  which  is  commonly  observed  during  the 
experimental  probing  of  the  NDR  current-region  of  RTD’s  is 
the  development  of  oscillations  with  bistable  behavior.  In  some 
cases,  these  spurious  oscillations  can  not  be  eliminated  and 
interfere  with  the  measurement  of  the  negative  conductance 
of  the  device  [12].  While  there  has  been  some  debate  over 
whether  these  oscillations  are  a  result  of  extrinsic  (i.e.  induced 
by  external  bias  circuit)  or  intrinsic  factors  [13],  [14],  recent 
quantum  mechanical  simulations  of  a  completely  isolated 
structure  indicate  the  RTD  has  an  inherent  ability  to  self- 
oscillate.  These  theoretical  results  [15]  are  free  of  the  extrinsic 
influences  always  present  during  experimental  measurements 
yet  exhibit  characteristics  observed  in  actual  measurements 
such  as  a  plateau  and  a  hysteresis  effect.  In  fact,  comparing 
the  I-V  results  from  these  quantum-based  simulations  (i.e. 
see  Fig.  2)  to  experimental  measurements  of  a  similar  RTD 
structure  (e.g.  [5])  reveals  a  remarkable  agreement  in  the 


1 


Stries-Inductance  Model  PanUel-Inductance  Model 

(a)  <*» 

Fig.  1.  Basic  circuit  configurations  used  in  nonlinear  oscillation  study. 


Applied  Bias  (V) 

Fig.  2.  Overall  current-density  versus  applied  bias  derived  from  Wigner 
distribution  simulation. 


plateau  and  step-like  features.  Hence,  these  observations  venfy 
the  presence  of  multiple  electric  storage  effects  in  the  intrinsic 
RTD  structure.  These  self-oscillations  are  useful  because  they 
can  be  used  as  a  guide  to  determining  an  accurate  equivalent 
circuit  model  for  the  purely  intrinsic  RTD. 

The  observation  of  self-oscillations  in  RTD’s  is  also  im¬ 
portant  from  another  perspective.  For  example,  previous  lim¬ 
itations  on  the  output  power  of  RTD-based  sources  can  be 
directly  linked  to  the  manner  of  their  implementation.  Tra  i- 
tionally  RTD’s  are  implemented  purely  as  a  NDR  element 
with  one  energy  storage  element  (i.e.  device  capacitance)^ 
Using  this  design  approach,  oscillauons  must  be  generated 
by  limit-cycles  which  exchange  energy  with  resonating  ele¬ 
ments  residing  in  the  external  bias  circuit.  This  approach  of 
extrinsically  inducing  oscillations  will  always  encounter  outpu 
power  restrictions  by  external  losses  (e.g.  contact  resistance) 
and  low  frequency  design  constraints  (i.e.  suppression  o 


bias  oscillations).  However,  this  theoretically-observed  self¬ 
oscillation  process  [15],  which  occurs  due  to  the  separate 
char^in 2-delav  times  of  the  emitter  accumulation  region  and 
the  quantum-well,  provides  evidence  for  a  possible  intrinsic 
approach  to  high-frequency  power  generation^  Since  these 
oscillations  have  frequencies  extending  up  to  6  THz  and  their 
amplitudes  are  not  limited  by  external  constraints,  there  is 
much  motivation  for  better  understanding  this  intrinsic  process. 
Hence,  the  identification  of  an  accurate  equivalent  circuit 
model  for  the  RTD  is  valuable  for  providing  insight  on  ways  to 
encourage  the  occurrence  of  this  intrinsic  phenomenon  and  for 
guiding  methods  to  couple  this  intrinsically-generated  power 
to  output  circuitry. 

This  paper  presents  a  comprehensive  evaluation  of  the 
circuit  model  configurations  presently  in  use  for  approximating 
the  behavior  of  RTD’s  and  extends  them  to  arrive  at  a  new 
model.  In  the  studies  presented,  the  numerically-generated 
nonlinear  results  associated  with  a  number  of  autonomous 
second-order  circuits  are  compared  to  the  self-oscillauons  and 
the  average  current-density  characteristics  observed  in  the 
quantum  mechanical  simulations  for  a  specific  case.  D^work 
identifies  a  specific  circuit  model  configuration  which  has 
the  potential  for  matching  the  current-density  characteristics 
observed  both  experimentally  and  theoretically.  Specifically, 
this  work  reveals  two  major  requirements  of  the  model  before 
it  can  accurately  reproduce  the  fundamental  self-oscillauons 
and  the  overall  I-V  characteristics  (i.e.  plateau  structure  and 
hysteresis).  First,  the  RTD  circuit  model  must  incorporate  a 
quantum-well  inductance  which  directly  chokes  the  nonlin¬ 
ear  conductance  element.  Second,  the  model  must  include 
a  nonlinear  “access”  resistance.  This  new  resistance  mode , 
associated  with  the  accumulation  of  charge  m  the  adjacen 
injection  (emitter)  region,  possesses  a  nonlocal  dependence 
on  the  bias  across  the  double  barrier  structure.  Die  results 
from  the  numerical  simulation  study,  presented  in  Section  IU, 
will  confirm  the  need  for  both  of  these  elements  in  the 
equivalent  circuit  model.  Specifically,  a  nonlinear  (i.e.  with 
nonlocal  voltage  dependence)  access  resistance  is  required  to 
match  the  known  hysteresis  behavior.  Furthermore,  the  circuit 
model  must  also  include  a  quantum-well  inductance  before 
it  can  produce  the  correct  self-oscillations  and  cunent-plateau 
phenomena.  Finally,  this  work  will  show  that  nonlinear  effects 
associated  with  this  inductive  element  are  responsible  for  some 
of  the  finer  features  of  the  oscillations  and  bistable  behavior. 
These  results  are  important  because  this  new  model  will  allow 
for  the  accurate  optimization  of  output  power  m  RTD-based 
high-frequency  sources  in  the  future. 


II,  Circuit  Model  Characteristics 
The  theoretical  identification  of  an  accurate  large- signal 
equivalent-circuit  model  for  the  resonant  tunneling  diode 
(RTD)  is  important  from  a  number  of  perspectives, 
accurate  circuit  model  is  important  in  oscillator  design  on 
(i)  maximizing  the  power-frequency  product  which  includes 
the  suppression  of  low-frequency  (LF)  oscillations  m  the 
bias  circuitry  [11],  [16];  (ii)  extracting  individual  component 
characteristics  from  experimental  measurements  [17]  !  J, 
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and  (iii)  providing  insight  for  structural  changes  to  the  RTD 
which  will  lead  to  an  improved  device-cavity  system  [6]. 

In  selecting  the  basic  RTD  circuit-model  configurations  for 
our  study,  we  will  utilize  its  fundamental  physical  charac¬ 
teristics  as  a  guide.  The  intrinsic  RTD  is  known  to  possess 
the  following  general  characteristics:  (i)  a  low-frequency  I-V 
characteristic  with  a  NDC  region  which  arises  from  tunneling 
through  quantum  states  near  the  resonant- state-energy  of  the 
double  barrier  structure:  (ii)  a  capacitance  due  to  the  negative 
charges  in  the  accumulated  emitter  and  well  imaging  of  the 
positive  charge  in  the  depleted  collector;  (iii)  an  access  resis¬ 
tance  associated  with  the  highly  doped  bulk  regions  on  either 
side  of  the  double  barrier  structure;  and  (iv)  a  self-inductance 
which  models  the  intrinsic  delay  in  the  tunneling  current 
with  respect  to  the  barrier-well-barrier  voltage  [15],  [20],  The 
basic  circuit  configurations  selected  for  this  nonlinear  study 
are  given  in  Fig.  1(a)  and  (b).  Since  every  fundamental¬ 
mode  oscillator  requires  at  least  two  energy-storage  elements 
[21],  these  two  autonomous  second-order  nonlinear  circuits 
represent  the  simplest  forms  which  simultaneously  satisfy  all 
the  characteristics  itemized  in  (i)  through  (iv).  The  series- 
inductance  (i.e.  inductance  is  in  series  with  device  capacitance) 
model  given  in  Fig.  1(a)  was  originally  proposed  by  Gering 
ex  al  [22]  to  explain  an  experimentally  observed  resonance 
matching  a  series-inductance  model.  The  parallel-inductance 
(i.e.  inductance  in  parallel  with  device  capacitance)  model 
given  in  Fig.  1(b)  is  identical  in  circuit-form  with  the  one 
proposed  by  Brown  et  al.  [20]. 

The  presentation  which  follows  will  compare  the  nonlinear 
large-signal  behavior  of  each  circuit  to  self-consistent  quantum 
mechanical  simulations  of  a  specific  RTD  structure.  In  this 
study,  each  circuit  is  evaluated  for  its  ability  to  match  self- 
consistent  Wigner-distribution-based  simulations  of  a  specific 
RTD  geometry  previously  presented  by  Jensen  and  Buot  [23]. 
This  RTD  is  a  GaAs-AlGaAs  based  structure  with  30  A  barri¬ 
ers  and  a  50  A  well.  The  AlGaAs  potential  barriers  are  0.3  eV 
and  the  device  temperature  is  77  K.  The  entire  simulation 
region,  including  2  x  1018  cm3  doped  buffer  regions,  is  550 
A.  In  this  study,  dissipation  effects  were  included  through  both 
internal  device  scattering  and  open  boundary  conditions.  Note, 
the  very  short  simulation  “box”  employed  does  not  influence 
the  results  since  “free-flowing”  (as  opposed  to  traditional 
ohmic)  boundary  conditions  were  enforced  [24],  [25].  Time- 
dependent  solutions  for  the  response  of  this  RTD  structure 
yield  the  overall  curTent-voltage  (I-V)  characteristic  displayed 
in  Fig.  2.  For  the  forward  sweep,  current  oscillations  were 
observed  within  the  voltage  region  V  =  0.24  to  V  =  0.32. 
For  the  backward  sweep,  current  oscillations  were  observed 
only  over  the  small  region  V  =  0.24  to  V  —  0.25.  Also, 
given  in  Fig.  2  is  a  plot  of  the  I-V  curve  when  the  entire 
applied  voltage  is  assumed  to  drop  linearly  across  the  quantum 
barrier-well  region.  For  this  linear  voltage-drop  model  all  the 
results  were  static  [15]. 

Note,  the  theoretical  I-V  curves  in  Fig.  2  display  peak 
current-densities  of  approximately  6  x  105  A/cm  2  which  are 
considerably  greater  than  that  typically  reported  for  RTD’s 
of  comparable  dimensions  in  this  material  system.  However, 
it  should  be  noted  that  conventional  GaAs-AlGaAs  RTD 


structures,  which  typically  possess  peak  densities  of  less  than 
1  x  105  A/cm“2.  almost  always  employ  undoped  spacer  layers 
(i.e.  between  buffer  and  barrier  regions)  and  usually  utilize 
limited  doping  levels  in  the  emitter  region  [12].  This  type  of 
design  is  chosen  for  practical  experimental  devices  to  insure 
undoped  barrier/well  regions  and  to  improve  device  contacts. 
Since  the  peak  current  density  of  an  RTD  is  highly  dependent 
on  the  electron  density  that  resides  at  the  emitter-barrier 
interface  [5],  these  relatively  higher  values  of  current-density 
obtained  for  the  structure  under  consideration  (i.e.  with  highly 
doped  emitter  and  no  spacer  layer)  should  be  expected. 

III.  Numerical  Simulation  Study 

In  this  section,  a  numerical  simulation  study  is  used  to  assess 
the  accuracy  of  circuit  models  presently  in  use  for  predicting 
the  high-frequency  performance  of  RTD’s.  Here,  a  van  der 
Pol  based  method  [26]  is  used  to  solve  for  the  nonlinear 
oscillations  associated  with  each  circuit  under  consideration. 
This  work  will  show  that  previous  models  lack  essential 
elements  necessary  for  predicting  the  self-oscillations  observed 
experimentally  during  measurements  of  the  current-voltage 
characteristics  [27]  and  predicted  theoretically  by  quantum 
mechanical  simulations  [15]. 

A.  Series-Inductance  Model:  Constant  Parameters 

The  first  equivalent  circuit  model  considered  is  the  series- 
inductance  model,  as  given  by  Fig.  1(a),  with  constant  values 
of  energy  storage  elements  (i.e.  L  and  C)  and  constant  access 
resistance  R .  Assuming  that  the  series- inductance  circuit  is 
biased  by  a  dc  voltage  source  Vs,  application  of  standard 
circuit  relations  readily  yield  the  autonomous  second-order 
nonlinear  equation, 

Vac)*^  +  v2{v,Vdc)v  =  0  (1) 

where 

i?(u,  Vdc)  =  {LG(v,Vd c)  4*  RC)/(LC) 

w2{v,  Vdc)  =  {Rg(v,Vd c)  +  1  )/(LC). 

Equation  (1)  results  from  a  transformation  into  the  dc  current- 
voltage  frame  of  reference  where  the  total  voltage  across 
the  nonlinear  element  is  V(t)  =  v(t)  4-  Vdc  and  VdC  is  the 
oscillation-free  voltage  in  the  circuit.  Hence,  (1)  describes  the 
evolution  of  v  which  is  the  transient  component  of  voltage. 
This  transformation  explains  the  dc  voltage  dependence  of  the 
parameters,  g(v,  Vdc)  =  {I(v  4-  Vdc)  ~~  I(Vd c))/v  which  can 
be  defined  as  the  nonlinear  conductance,  since  in  this  frame 
the  transient  current  density  is  i  =  I{v  4-  Vdc)  —  I(Vdc)  = 
g(v ,  Vdc)v,  and  G(v,  Vdc)  =  G{ V)  =  dl/dV  which  is  just  the 
slope  of  I(V)  at  V  =  v  4-  Vdc  or  the  normal  definition  of  the 
linearized  conductance.  This  transformation  allows  solution  of 
the  nonlinear  dc  biasing  separately  from  the  time-dependent 
oscillation  problem,  although  dynamic  biasing  effects  are  still 
inherently  included. 

After  a  modified  Newton’s  method  [28]  was  used  to  solve 
the  biasing  problem,  a  numerical  version  of  the  van  der  Pol 
algorithm  [29]  was  employed  to  solve  for  the  unstable  (or 
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Fig.  3.  Overall  current-density  solutions  for  forward,  backward  and  tri-stable 
state  for  condition  of  absolute  stability  (i.e.  no  self-oscillations). 

decaying)  states  of  (1).  This  van  der  Pol-based  algorithm  is 
simply  two  successive  numerical  integrations  performed  on  the 
system  which  results  when  (1)  is  written  as  two  first  order 
equations  [26].  Specifically,  introduction  of  the  substitution 
variable  2  =  dv/dt  =  v  and  using  2  =  z  ■  ( dz/dv )  leads  to 
the  trajectory  equation 

fli  =  -ri(v,Vdc)  -  w(«,Vdc)(-)-  (2) 

In  analogy  with  the  method  of  isoclines  or  the  phase-portrait 
technique  [21],  (2)  is  numerically  integrated  to  generate  a 
phase-space  trajectory  for  2  as  a  function  of  v.  Solution  curves 
for  v  as  a  function  of  time  are  then  obtained  by  approximate 
integration  of  2  =  dv/dt. 

For  the  initial  circuit  simulations  of  this  case  study,  param¬ 
eter  values  were  estimated  from  previous  quantum  mechanical 
calculations.  The  nonlinear  conductance  I(V)  is  defined  from 
the  linear  voltage-drop  model  result  in  Fig.  2.  The  constant 
value  for  R  =  1.63  x  10“7  ft  cm2  was  approximated  from 
the  shift  in  the  peak  of  the  I-V  characteristic  given  in  Fig.  2. 
The  constant  value  of  C  =  2.3  x  10-7  F/cm2  was  obtained 
by  calculating  the  total  displaced  charge  as  a  function  of  the 
bias  voltage  and  taking. the  average  over  the  NDR  region.  The 
value  of  L  =  1.48  x  10"21  H  cm2  was  estimated  by  matching 
a  linear  oscillation  analysis  to  the  current-oscillations  observed 
during  the  physical  simulations.  For  a  detailed  discussion  of 
the  reasoning  for  these  initial  parameter  choices  refer  to  [15]. 

The  first  step  of  the  analysis  is  the  determination  of  V^c 
which  is  obtained  from  the  solution  of  I{V)  +  ( V- Vs)/ R  =  0 
where  Vs  is  the  applied  bias.  The  important  point  to  note 
is  that  for  some  values  of  Vs  there  are  multiple  solutions 
for  the  Vic -intercept.  For  example,  consider  the  complete  dc 
mapping  which  results  from  all  intercept  points  as  given  in 
Fig.  3.  Clearly,  a  hysteresis  is  present  in  this  oscillation-free 
current- voltage  characteristic  Id(Vs)  where:  (i)  the  forward 
voltage  sweep  solution  is  associated  with  the  low- Vs  intercept 
point;  (ii)  the  backward  voltage-sweep  solution  is  associated 
with  the  high- V5  intercept  point:  and  (iii)  the  tri-stable  solution 
is  associated  from  the  midpoint  intercept. 


The  van  der  Pol  algorithm  was  then  applied  to  both  the 
forward  and  backward  sweeps  to  study  the  general  stability 
of  the  series-inductance  circuit  model.  The  introduction  of 
noise- vectors  (impulses)  and  the  subsequent  generation  of 
the  z-v  trajectories  showed  decaying  solutions  for  ail  Vdc 
intercepts  both  within  and  outside  the  NDR  (i.e.  region  where 
G(v,  Vdc)  <  0)-  The  same  procedure  was  applied  to  the  tri¬ 
stable  VdC  intercepts  and  in  all  cases  transients  were  observed 
which  either  converged  to  the  forward-sweep  or  the  negative- 
sweep  Vdc  intercepts  appropriate  for  the  applied  bias  under 
consideration.  A  variational  study  was  also  performed  to 
determine  if  values  of  L  and  C  could  be  found  which  led  to 
oscillations.  Studies  over  a  wide  range  (i.e.  ten  times  larger  and 
smaller)  of  the  base  values  Co  =  2.3  x  10“ '  F/cm2  and  Lo  — 
1.48  x  10~21  H  cm2  yielded  the  same  stable  results.  These 
results  agree  with  a  previous  linear  analysis  which  indicated 
that  the  series-inductance  model  was  incapable  of  yielding 
the  self-oscillations  observed  in  the  physical  simulations  with 
these  estimated  circuit  parameters  [15]. 

In  an  additional  variational  study  of  the  access  resistance 
R ,  oscillations  were  found  to  exist  for  values  R  <  0.5  x 
10~7  ft  cm2.  Fig.  4  shows  the  phase-space  trajectory  of  the 
main-lead  current  density  for  R  =  0.5  x  10“7  ft  cm2  and  L  — 
Lo  =  1.48  x  10“21  H  cm2  at  an  applied  bias  of  V5  =  0.179  V. 
The  period  of  this  oscillation,  which  is  0.15  ps,  matches  one 
of  the  self-oscillations  of  the  quantum  simulations;  however, 
it  occurs  at  a  much  lower  applied  bias.  Fig.  5  compares  the 
average  forward-sweep  current  versus  bias  for  the  three  cases 
L  =  Lo/5Lo  md  5Lo  with  R  =  0.5  x  1<T7  ft  cm2.  For  the 
oscillating  cases,  there  is  a  development  of  a  current-plateau 
for  increasing  values  of  L  and  a  similar  result  (not  shown)  can 
be  obtained  for  variations  in  C.  For  clarity,  we  have  not  shown 
the  backward-sweep  results  which  exhibits  a  hysteresis  similar 
in  extent  to  that  seen  in  Fig.  3.  These  general  results  agree  with 
previous  studies  of  circuits  with  the  series-inductance  form 
[30],  [31  ].  The  main  point  to  note  is  while  one  can  achieve  self¬ 
oscillations  and  a  hysteresis  effect  with  the  constant-parameter 
series-inductance  model  by  properly  choosing  L  (or  C)  this 
can  only  be  done  by  significantly  reducing  the  estimated  value 
of  access  resistance  R.  Since  this  reduction  of  R  leads  a  large 
backward  shift  in  the  overall  Id{Vs)  characteristic,  an  overall 
match  to  the  physics-based  simulations  is  not  possible. 

B.  Parallel-Inductance  Model:  Constant  Parameters 

The  next  circuit  model  considered  is  the  parallel-inductance 
circuit  model  with  constant  parameters  as  given  in  Fig.  1(b). 
For  the  investigations  of  this  model,  the  initial  values  em¬ 
ployed  for  R.  L  and  C  will  be  the  same  as  those  defined 
in  the  last  section.  The  dynamics  of  the  parallel-inductance 
circuit  are  described  by  the  nonlinear  equation 

$+^,vde)|+C(«Tdc)(|)2W(t,,vdc)V  =  0(3) 

where 

Kic)  —  (LG(v,  Vdc)  +  RC)/(RCLG(v.  Vdc)) 
C(«!Vdc)  =  (G>.Vdc))/(G(v.Kic)) 
tjj2(v,  Vdc)  =  (1  +  Rg(v,  VdC))/(RCLG(v,Vd  c)). 
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Fig.  4.  Current-voltage  phase-space  trajectory  of  a  self-oscillating  condi¬ 
tion 


Fig.  5.  Overall  current-density  versus  applied  bias  derived  from  se¬ 
ries-inductance  circuit  model  with  constant  parameters. 


Here,  as  before,  g(v,  Vdc)  is  the  nonlinear  conductance  and 
G{v.Vdc)  is  the  linear  conductance.  Since  the  dc  circuit  for 
this  model  is  the  same  as  that  for  the  series-inductance  model, 
the  Vac  reference  is  defined  by  the  same  loadline  intercepts 
as  in  the  previous  section.  Also,  note  the  emergence  of  the 
new  parameter  C^^dc)  containing  the  relation  G'(v,  Vdc)  = 
G'(V)  =  dG/dV  which  is  just  the  derivative  of  the  linear 
conductance.  Unlike  the  defining  equation  for  the  previous 
model,  the  coefficients  for  (9)  all  contain  G{v,Vic)  in  their 
denominators;  and  hence,  possess  the  potential  to  become 
singular.  Later,  this  will  be  shown  to  lead  to  “jumps”  in  the 
phase-space  tniectories. 

The  numerical  algorithm  was  employed  to  investigate  the 
stability  of  the  parallel-inductance  model  for  forward  and 
backward  bias  sweeps.  After  determining  the  appropriate  Vdc- 
intercepts,  the  z-v  space  trajectory  is  then  obtained  by  inte¬ 
gration  of 

=  -T,(v,  Vdc)  -  <(t).  vdc)z  -  W2(t).  Vdc)Q  (4) 


Fig.  6.  Overall  current-density  versus  applied  bias  derived  from  paral¬ 
lel-inductance  circuit  model  with  constant  parameters. 


and  a  time  mapping  is  then  achieved  from  z  =  dv/dt.  Fig.  6 
gives  the  overall  current-density  results  for  the  estimated 
circuit  parameters  R  —  1.63  x  10  1  $7  cm2.  C  =  2.3  x 
10~7  F/cm2,  and  L  =  1.48  x  10“21  Hem2.  Note,  the  dc 
current  components  are  also  included  for  reference.  In  the 
applied  bias  range  (0.20.0.23)  the  characteristic  is  found 
to  display  a  distinct  hysteresis  with  oscillations  in  both  the 
forward  and  backward  sweeps  within  the  NDR.  The  origin 
of  this  hysteresis  can  be  understood  by  referring  to  Fig.  7 
which  shows  the  phase-space  trajectories  for  forward  and 
backward  bias  sweeps  at  Vs  =  0.22  V.  While  the  final  limit- 
cycles  of  the  two  trajectories  are  equivalent,  the  dc  bias  points 
for  the  forward  and  backward  differ  by  a  significant  value. 
Thus,  the  overall  (or  total)  current-density  is  a  function  of 
sweep  direction.  Another  feature  of  note  is  the  very  rapid  (but 
finite)  changes  in  the  current-density  at  the  phase-space  turning 
points.  Fig.  8  gives  the  time  domain  representation  of  the 
forward-sweep  trajectory  from  Fig.  7.  As  shown,  these  jumps 
in  current-density  correspond  to  instances  where  dv/dt  is 
discontinuous.  These  allowable  discontinuities  in  dv/dt  arise 
from  the  intersection  of  the  trajectory  with  the  singular  points 
defined  by  G(v,  Vdc)  =  0  (recall  this  drives  the  coefficients 
in  (9)  to  infinity).  This  phenomenon  can  be  understood  by 
referring  to  Fig.  9  which  shows  the  final  z-v  space  limit-cycle 
for  the  forward-sweep  trajectory  from  Figs.  7  and  8.  As  the 
value  of  v  approaches  the  point  where  G(v .  Fdc)  =  0,  the  z 
ordinate  is  seen  to  approach  infinity.  However,  a  limit  analysis 
of  (9)  as  G  —  0  reveals  that  it  degenerates  to  the  first-order 
quadratic  equation 


(z)2  + 


(  1  ,  ( Rg(vj,v dc)  +  iy 

\G’(vj.Vdc)  )Z  +  \RCLG(vJ:Vdc)  J 


=  0 


(5) 


where  z  =  dv/dt  and  vj  is  the  jump  value.  Hence,  as  the 
z-v  space  trajectory  makes  its  clockwise  cycle  and  these 
singularities  are  encountered  a  jump  occurs  at  the  turning 
points.  Here,  the  jump  is  from  z  —  oc  to  the  appropriate 
root  of  (5).  In  the  numerical  algorithm,  all  that  is  required  is 
to  estimate  the  location  of  the  v-space  turning  point(s)  and 
to  integrate  to  a  sufficiently  large  value  of  z  to  ensure  the 
accuracy  of  the  v  and  t  mapping. 
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Fig.  7.  Current-density  phase-space  trajectories  for  forward  and  backward 
bias  sweeps  at  an  applied  bias  of  0.22  V. 


Time  (ps) 

Fig.  8.  Time  evolution  of  the  total  current-density  for  the  forward  bias  sweep 
at  an  applied  bias  of  0.22  V.  Here,  the  internal  conductance  bias  (scaled  by 
a  factor  of  ten  is  volts)  is  also  plotted  for  reference. 

Returning  to  a  discussion  of  Fig.  6,  the  hysteresis  collapses 
at  Vs  ~  0.23  V;  however,  a  plateau  remains  out  to  Vs  =  0.27 
V.  This  plateau  is  the  result  of  an  oscillation  about  an  attractor 
in  z-v  space  which  occurs  at  a  point  near  the  low-voltage 
singularity.  The  oscillations  of  this  plateau  region  are  of  a 
very  high  frequency  (period  <  0.01  ps)  and  of  very  small 
amplitude  (v  <  0.01  x  105  A/cm2).  Thus,  these  current- 
density  characteristics  would  appear  as  static  if  they  were 
derived  from  simulations  in  a  total  reference  frame  or  from 
actual  measurements. 

Simulations  were  also  performed  to  determine  how  vari¬ 
ations  in  L  anchor  C  would  influence  the  overall  results. 
While  these  studies  indicate  that  the  hysteresis  and  plateau 
structure  can  be  significantly  altered  by  changes  to  these 
terms,  the  major  discrepancies  between  the  results  of  Fig.  6 
and  the  physical-simulation  results  of  Fig.  2  could  not  be 
resolved.  Limitations  of  the  parallel-inductance  model  with 
constant  parameters:  (i)  a  restricted  region  of  hysteresis;  (ii) 
oscillations  present  in  the  backward  sweep;  and  (iii)  incorrect 
amplitude  and/or  frequency  of  forward  sweep  self-oscillations. 


Fig.  9.  Limit-cycle  of  the  z-v  space  trajectory  for  the  forward  bias  sweep  at 
an  applied  bias  of  0.22  V.  This  limit-cycle  exhibits  jumps  at  the  turning  points 
which  correspond  to  the  peak  and  valley  of  the  nonlinear  I(  V)  characteristic. 


Therefore,  the  constant-parameter  parallel-inductance  model  is 
not  suitable  for  describing  the  dynamics  of  the  specific  RTD 
under  consideration. 

The  results  of  the  last  two  sections  have  shown  that  neither 
of  the  basic  circuit  configurations  with  constant  parameters 
is  sufficient  to  describe  the  observed  self-oscillations.  While 
one  should  expect  the  presence  of  nonlinearities  in  the  storage 
elements  (i.e.  C  and  L),  the  previous  results  suggest  that  the 
constant  R  is  a  barrier  to  achieving  the  correct  hysteresis 
in  the  overall  current-density  characteristic.  The  value  of  the 
constant  R ,  used  up  to  this  point,  was  estimated  from  the  shift 
in  the  RTD’s  I-V  characteristic  when  self-consistency  effects 
are  included.  This  approach  is  accurate  only  if  the  true  access 
resistance  is  independent  of  bias.  Since  the  backward  sweep 
current-density  given  in  Fig.  2  is  almost  (i.e.  except  for  the 
small  bias  region  (0.24.0.25))  free  of  oscillations  the  voltage 
dependency  of  R  can  be  easily  estimated.  A  point  by  point 
estimation  of  the  shift  in  the  current-density  can  be  used  to 
determine  the  effective  resistance  as  a  function  voltage  across 
the  barrier  structure.  The  character  of  this  voltage -dependent 
resistance  R(V)  is  given  approximately  by  the  piece-wise 
curve  shown  in  Fig.  10.  The  influence  of  R{V)  on  the  previous 
circuit  models  will  be  considered  in  the  sections  which  follow. 

C  Series- Inductance  Model :  Nonlinear  R(V)  Model 

The  transient  response  of  the  series -inductance  model  with 
nonlinear  access  resistance  R{V)  and  constant  storage  ele¬ 
ments  is  described  by  the  nonlinear  equation 

<^+v(v,Vdc)^+u2{v,Vdc)v  =  0  (6) 

where  the  coefficients  are  the  same  as  in  (1)  with  R{V) 
substituted  for  R.  Here,  as  before,  V^c  is  the  reference  frame 
which  is  found  now  by  solving  the  dc  loadline  equation 
with  the  nonlinear  J?(V).  Fig.  11  displays  a  number  of  the 
loadline  intercepts  for  (15)  where  R(V)  is  defined  from  the 
characteristic  given  in  Fig.  10.  This  nonlinear  access  resistance 
causes  the  resultant  loadlines  to  sag  downward  in  the  NDR 
region  increasing  the  range  of  Vs  values  over  which  multiple 
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Fig.  10.  Access  resistance  as  a  function  of  voltage  across  the  intrinsic  RTD 
structure. 


Fig.  11.  Intersection  curves  of  I{V)  characteristic  and  nonlinear  loadiine 
for  a  range  of  applied  biases. 


V^-intercepts  occur.  Fig.  12  plots  the  overall  current-density 
characteristic,  for  forward  and  backward  bias  sweeps,  obtained 
from  using  the  base  circuit  parameters  Co  =  2.3  x  10” '  F/cm2 
and  L0  =  1-48  x  1CT21  Hem2.  In  the  forward  sweep,  oscilla¬ 
tions  are  obtained  in  the  voltage  range  (0.23.0.30);  however, 
these  oscillations  led  to  a  very  large  downward  shift  in 
the  average  current.  In  the  backward  sweep,  the  circuit  was 
stable  over  most  of  the  hysteresis  region  with  an  up-shifting 
oscillation  occurring  over  the  voltage  range  (0.23,0.25).  For 
both  sweep  directions  the  oscillations  result  from  a  process 
of  cycling  closely  about  a  phase-space  attractor  point;  hence, 
they  are  of  a  very  high  frequency  (period  <  0.01  ps)  and 
small  amplitude  (v  <  0.01  x  105  A/cm2).  In  addition,  the 
average  current  results  were  found  to  be  entirely  independent 
of  the  values  of  L  and  C.  This  result  differs  greatly  from 
a  purely  linear  analysis  [15]  which  predicts  stable  results  if 
RC  <  L\G\.  Since  it  is  not  possible  to  match  either  the 
forward-sweep  plateau  height  or  the  nature  of  the  oscillations 
observed  in  the  physics-based  simulations  this  model  is  an 
insufficient  description  of  the  RTD. 


Fig.  12.  Overall  current-density  versus  applied  bias  derived  from  sc- 
ries-inductance  circuit  with  nonlinear  R{V)  model. 


D.  Parallel-Inductance  Model:  Nonlinear  R{V)  Model 

The  transient  response  of  the  parallel-inductance  model 
with  nonlinear  access  resistance  R(V)  and  constant  storage 
elements  is  described  by  the  nonlinear  equation 

v^Tt +av- Vdc){jt )2  +oj2(v'Vi‘)v  =  0  (7) 

where  the  coefficients  are  the  same  of  those  of  (3)  with 
R  replaced  by  R{V ).  Here,  as  in  the  last  section,  V<±c  is 
found  by  solving  the  nonlinear  loadline  equation  (i.e.  see 
Fig.  12).  The  van  der  Pol  algorithm  was  applied  to  this 
form  of  the  parallel-inductance  model  to  determine  the  phase- 
space  trajectories.  Fig.  13  displays  the  overall  current-density 
results,  for  forward  and  backward  sweeps,  obtained  from  using 
the  base  circuit  parameters  Co  =  2.3  x  10” '  F/cm2  and 
Lo  =  1.48  x  10“21  Hem2.  These  results  indicate  a  very  good 
match  between  the  characteristic  obtained  from  the  circuit 
model  and  those  of  the  physical  simulations.  For  example,  a 
forward-sweep  plateau  is  present  with  correct  magnitude  (i.e. 
«4.5  x  105  A/cm2)  and  extent  (  i.e.  V  =  0.23-0.30).  Though 
not  explicitly  apparent  in  the  average  results  of  Fig.  13, 
the  correspondence  between  the  presence  of  oscillations  and 
applied  bias  almost  exactly  matched  the  physical  simulations 
in  both  the  forward  and  backward  sweep.  Unfortunately,  for 
the  base  inductor  value  Lo  the  oscillations  exhibited  very 
large  peak-to-peak  swing  and  rapid  jumps  in  the  current- 
density  at  the  phase-space  turning  points.  However,  when  the 
inductor  value  was  changed  to  L  =  5.5  x  10”21  Hem2  to 
mne  the  oscillation  amplitude  to  the  physics-based  results  at 
Vs  =  0.26  V,  a  very  good  frequency-match  was  achieved 
over  the  entire  plateau-region.  For  example.  Fig.  14  gives  the 
time  evolution  of  the  current-density  at  two  different  bias 
using  this  new  value  of  L.  As  shown,  the  periods  are  now 
approximately  0.3  ps  and  0.17  ps  at  Vs  =  0.26  V  and  0.28  V 
compared  to  quantum  simulation  results  of  0.4  ps  and  0.133  ps, 
respectively.  While  the  dependance  of  oscillation  amplitude  on 
bias  was  not  an  exact  match,  these  results  identify  the  parallel- 
inductance  circuit  model,  with  nonlinear  R{V)  element,  as  an 
excellent  candidate  for  describing  the  dynamics  of  the  RTD. 
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Fig.  13.  Overall  current-density  versus  applied  bias  derived  from  paral¬ 
lel-inductance  circuit  with  nonlinear  /2(V')  model. 


Fig.  14.  Time  evolution  of  the  total  currem-density  for  the  forward  bias 
sweep  at  an  applied  bias  of  0.26  and  0.28  V  with  L  =  5.5  x  10-21  Hem2. 
Here,  the  internal  conductance  biases  (scaled  by  a  factor  of  ten  in  volts)  are 
also  plotted  for  reference.  Note,  upper  curves  are  current-densities  and  lower 
curves  are  internal  bias. 

E.  Parallel-Inductance  Model:  Nonlinear 
R(V)  and  L{V)  Models 

The  previous  studies  have  shown  that  the  parallel- 
inductance  model  with  a  nonlinear  access  can  successfully 
match  many  of  the  characteristics  associated  with  the  self¬ 
oscillations  of  RTD’s.  The  only  real  weakness  of  the  model 
is  its  inability  to  simultaneously  tune  to  the  frequency  and 
amplitude  variations  as  a  function  of  applied  bias.  Since  the 
value  of  L  directly  modifies  the  growth/decay  coefficients  of 
the  nonlinear  transient  equation  given  in  (7),  the  introduction 
of  a  nonlinear  inductance  seems  to  be  one  possible  approach 
of  overcoming  the  limitations  of  the  circuit  model  considered 
in  the  last  section.  In  fact  Brown  et  al  [20]  has  previously 
proposed  the  introduction  of  a  quantum  well  inductance  given 
by  Lq w  —  r/G  where  r  is  the  quasibound  state  lifetime  and 
G  is  the  linear  conductance,  to  describe  the  delay  associated 
with  resonant  tunneling.  Since  this  nonlinear  study  utilizes 
a  nonlinear  I(V)  characteristic,  consistency  requires  the 
inductor  model  be  of  the  form  r/G(V)  which  in  our  dc 
reference  frame  becomes  r/G(v.  Vdc).  Note,  when  this  model 


Fig.  15.  Overall  current-density  versus  applied  bias  derived  from  paral¬ 
lel-inductance  circuit  with  nonlinear  R(Y)  model  and  nonlinear  L(V)  model. 


is  utilized  the  inductor  has  a  negative  value  for  points  inside 
the  NDR  region.  With  this  inductance  model,  the  dynamics 
are  described  by  the  nonlinear  equation 

»2  j 

+  v(v.  Vdc)^  +  u2(v.  Vdc)v  =  0  (8) 

v(v.Vdc)  =  l/r  +  l/(R(v.Vdc)C) 
where 

u2(v,  Vtc)  =  (R(V,  Vdc)g(v ?  vdc)  +  1  )/(tR(v.  Vdc)C). 

Notice  that  the  introduction  of  this  L(V)  model  actually  results 
in  a  simplified  version  of  (7)  where  now  the  new  coefficients 
in  (8)  are  not  able  to  become  singular.  The  initial  studies 
utilized  the  base  circuit  parameter  Co  —  2.3  x  10"7  F/cm2 
and  the  nonlinear  access  resistance  R{V)  defined  from  Fig.  10. 
In  addition.  The  initial  value  of  r  was  obtained  by  equating 
the  nonlinear  model  for  inductance,  at  the  midpoint  of  the 
NDR  region,  to  the  base  value  of  the  inductance  Lq  =  1.48  x 
1(T21  Hem2.  Utilizing  this  approach  yields  r  =  L0G(V  = 
0.15  V)  =  0.015  ps  which  is  in  approximate  agreement  which 
other  estimations  of  tunneling  time  associated  with  the  type  of 
RTD  under  consideration  [32].  Interestingly,  this  circuit  model 
yielded  results  identical  to  the  those  obtained  for  the  series- 
inductance  circuit  with  a  nonlinear  R(V )  model  studied  in 
Section  III-C.  Hence,  the  introduction  of  Brown's  inductance 
model  removes  the  singularities  that  bound  the  oscillation  to 
the  NDR  and  allow  the  trajectory  to  converges  to  a  phase- 
space  attractor  point.  As  in  Section  III-C.  the  overall  current 
results  of  this  circuit  were  found  to  be  independent  of  L  (i.e. 
r)  and  C. 

To  gather  additional  information  on  the  influence  of  nonlin¬ 
ear  inductance  on  the  parallel-inductance  model,  simulations 
were  also  performed  for  an  inductance  model  defined  from 
the  negative  of  Brown's  inductor  model.  Note,  this  results  in 
a  positive  inductance  inside  the  NDR  region  and  a  negative 
inductance  outside  the  NDR  region.  For  the  initial  studies 
utilizing  the  previous  base  parameters,  trajectories  for  both  the 
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forward  and  backward  were  found  to  become  unbounded.  To 
remedy  this  the  state  lifetime  was  increased  to  r  =  0.062  ps 
(through  a  variational  study)  which  yielded  bounded  limit 
cycles  for  the  forward  bias  sweep.  Fig.  15  shows  the  overall- 
current  density  for  the  parallel-inductance  model  with  nonlin¬ 
ear  access  resistance  and  negative  quantum-well  inductance 
model.  These  results  indicate  that  a  nonlinear  inductive  effect 
can  lift  the  higher-bias  portion  of  the  forward-sweep  current- 
density  plateau  which  is  a  characteristic  obser/ed  in  the 
physics-based  results  given  in  Fig.  2.  The  introduction  of  this 
horn  in  the  1{V)  curve  results  from  the  nonlinear  inductance 
which  admits  larger  values  of  L  for  positive  excursions  about 
the  Vdc- intercept  point.  Unfortunately,  in  the  unstable  part  of 
the  backward  bias  sweep  the  oscillations  were  still  found  to 
become  unbounded.  Also,  the  amplitude  of  the  forward-sweep 
oscillations  still  increases  with  increasing  bias.  However,  the 
results  of  this  section  have  shown  that  a  nonlinear  inductance 
element  can  be  used  to  improve  the  match  between  the 
equivalent  circuit  model  and  the  physics-based  results. 

IV.  Conclusion 

Numerically  generated  results  have  been  used  to  judge 
the  feasibility  of  a  number  of  equivalent  circuits  presently 
in  use  to  model  the  dynamics  of  resonant  tunneling  diodes 
(RTD’s).  In  these  studies  two  main  circuit  configurations  were 
considered,  each  of  which  possesses  the  general  electrical 
characteristics  of  RTD’s.  These  studies  have  verified  that 
neither  configuration  with  constant  R .  L,  and  C  elements  is 
sufficient  for  describing  the  self-oscillations  and  overall  current 
characteristics  observed  experimentally  and  predicted  theoreti¬ 
cally  by  fully  self-consistent  quantum  mechanical  simulations. 
These  studies  have  shown  that  the  introduction  of  nonlinear 
access  resistance,  with  nonlocal  dependence  on  the  bias  across 
the  NDR  current  source  (or  equivalently  the  intrinsic  double 
barrier  structure),  can  improve  the  match  to  the  hysteresis 
phenomena,  in  both  circuit  model  configurations.  However, 
only  the  circuit  with  a  quantum-well  inductance  which  directly 
chokes  the  nonlinear  conductance  element  is  capable  of  prop¬ 
erly  describing  the  current-plateau  phenomena.  In  addition, 
this  work  has  shown  that  incorporating  nonlinearity  effects  into 
this  inductive  element  can  resolve  some  of  the  finer  features  of 
the  oscillations  and  bistable  behavior.  While  additional  study  is 
required  to  derive  globally-accurate  expressions  for  describing 
the  full  nonlinearity  of  the  storage  elements  (i.e.  L  and  C\  this 
work  has  provided  an  improved  tool  for  the  optimization  of 
RTD-based  high  frequency  sources. 
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Group  on  Electron  Devices  (AGED)  and  the  DoD  JDL/Reliancc  Sub  Panel 
on  Solid  State  RF  Technology.  He  is  a  member  of  the  Army’s  Electronic 
Coordinating  Committee  (ECOG),  the  Technical  Advisory  Committee  (TAC) 
for  the  URI/High  Frequency  Microelectronics  Center  at  the  University  of 
Michigan  and  the  NSF/CAEME  policy  board  at  the  University  of  Utah. 
He  holds  four  U.S.  patents  and  has  published  over  65  technical  papers.  He 
has  received  several  awards  including  four  RCA  Outstanding  Engineering 
Achievement  Awards  for  research  in  microwave/millimeter-wave  device  and 
circuit  technology. 


David  L.  Rhodes  (S’80-M’80)  received  the  B.S.  degree  in  electrical  engi¬ 
neering  from  Rutgers  University,  New  Brunswick,  NJ,  in  1980  and  the  M.S. 
degree  from  Princeton  University.  Princeton,  NJ,  in  1982  while  working  at 
RCA  Laboratories,  David  Samoff  Research  Center,  Princton. 

He  joined  RCA  Laboratories  in  June  1980,  working  in  the  Microwave 
division.  There  he  primarily  worked  on  CAT  and  CAD  for  microwave 
applications  including  the  areas  of  test  automation,  simulation,  optimization, 
data  visualization,  and  statistical  design  techniques.  He  is  now  Branch  Chief  of 
the  Electrophysics  and  Modeling  branch  at  the  Physical  Sciences  Directorate 
(PSD)  of  the  Armv  Research  Laboratory  (ARL),  Fort  Monmouth.  NJ. 

He  is  chair  of  the  “IEEE  SCC-30  MHDL  Study  Group"  and  is  working  on  or 
directing  efforts  in  the  areas  of  parallel  simulation,  meta-modeling/statistical 

design,  and  active  device  modeling  in  addition  to  analog  HDL-related  work.  * 

He  is  currently  tn-service  leader  of  the  Computational  Electronics  and 
Nanoelectronics  area  of  High  Performance  Computing  Modernization  Program 
software  initiative.  He  was  the  technical  program  chair  for  the  first  and  second 

( 1994  &  1 995)  *' International  Conference  on  Electronic  Hardware  Description  * 

Languages  ( ICEHDLf  (called  SHDL  in  1994)  conferences.  Mr.  Rhodes  is 

a  member  of  the  IEEE  Circuits  &  Systems  and  Microwave  Theory  and 

Techniques  societies  and  the  Analog  VHDL  International  s  (AVI)  Executive 

Committee.  He  also  serves  on  the  Editorial  Board  for  the  International  Journal 

of  Microwave  and  Millimeter-Wave  Computer-Aided  Engineering. 
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